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The transport equations for polarized radiation transfer in non-Riemannian, Weyl-Cartan type 
space-times are derived, with the effects of both torsion and non-metricity included. To obtain the 
basic propagation equations we use the tangent bundle approach. The equations describing the time 
evolution of the Stokes parameters, of the photon distribution function and of the total polarization 
degree can be formulated as a system of coupled first order partial differential equations. As an 
application of our results we consider the propagation of the cosmological gamma ray bursts in 
spatially homogeneous and isotropic spaces with torsion and non-metricity. For this case the exact 
general solution of the equation for the polarization degree is obtained, with the effects of the torsion 
and non-metricity included. The presence of a non-Riemannian geometrical background in which 
the electromagnetic fields couple to torsion and/or non-metricity affect the polarization of photon 
beams. Consequently, we suggest that the observed polarization of prompt cosmological gamma ray 
bursts and of their optical afterglows may have a propagation effect component, due to a torsion/non- 
metricity induced birefringence of the vacuum. A cosmological redshift and frequency dependence 
of the polarization degree of gamma ray bursts also follows from the model, thus providing a clear 
observational signature of the torsional/non-metric effects. On the other hand, observations of the 
polarization of the gamma ray bursts can impose strong constraints on the torsion and non-metricity 
and discriminate between different theoretical models. 

PACS numbers: 95.30.Gv, 95.30.Sf, 98.70.Rz 



I. INTRODUCTION 



The study of the propagation of electromagnetic radiation in gravitational fields plays an essential role in relativistic 
astrophysics and cosmology. In order to obtain a correct and consistent description of the radiative transfer processes 
for photons propagating on cosmological distances or on the cosmic microwave background one must also take into 
account the effects of the background geometry on the radiation emitted near the surface of neutrons stars or black 
holes The radiation follows curved paths according to the underlying geometry and is red or blue-shifted. Spatial 
curvature in some geometries, like, for example, the Kerr geometry, rotates the polarization of photons propagating 
through them. 

Generally, one cannot solve Maxwell's equations exactly for waves propagating in curved space-times or in a relativis- 
tic medium. Instead one either uses a geometrical optics approximation or a kinetic description of the multi-photon 
system, by treating photons as massless classical particles, characterized by a four- momentum p at an event x 0]. 
In the geometrical optics approximation, which is reasonable in most of the cases, the degree of polarization is not 
affected by the underlying geometry. The field vector propagate parallel along the light paths. For non-flat curved 
paths, a rotation of the polarization angle results. 

Astrophysicists measure intensities and spectral distributions. These can conveniently be described by a distribution 
function, which is defined on phase space and is directly related to the spectral intensity. In this picture the propagation 
of radiation is described by an equation of radiative transfer, which is a differential equation for the photon distribution 
function. In addition to the above-mentioned quantities, degrees and plane of polarization of the observed radiation 
provide important means of information about the source and the intervening medium and geometry. 

A general relativistic form of the Boltzmann transport equation for particles or radiation interacting with an 
external medium has been developed by Lindquist 3] . The comoving frame radiation transfer equation in an Eulerian 



'Electronic address: hrspksc@hkucc.hku.hk 
^Electronic address: harko@hkucc.hku.hk 
t Electronic address: xywang@nju.edu.cn 



2 



coordinate system was obtained by Riffert |4|. A transport equation for photons in curved Riemannian space-times 
was derived by Bildhauer by making use of the covariant generalization of the Wigner transformation. Two linear 
equations on the tangent bundle are obtained whose correction terms to the Liouville equation and the classical mass- 
shell condition are of the order of the wavelength over the variation length of the background geometry. The comoving 
frame transfer equation for the Stokes parameters in an arbitrary background metric in an Eulcrian coordinate system 
has been obtained in Q and applied to the study of the Berry's phase for polarized light in Q. 

A long time ago Brans has suggested that a rotation of the polarization plane of the radiation arises whenever 
there is shear, by an effect analogous to the Thomas precession. Therefore if there is a large scale anisotropy in the 
expansion of the Universe, the microwave background radiation is expected to be linearly polarized. In Friedmann- 
Robertson- Walker models with small expansion anisotropy the observed rotation of the polarization plane would be 
appreciable and constant over the celestial sphere in closed (Bianchi type IX) models or it should vanish in flat (type 
I) and open (type V) models Q . Hence the study of the polarization of the cosmic microwave background radiation 
can imposes strong limits on the anisotropy and other fundamental physical parameters of the Universe. A rigorous 
first order solution with respect to the anisotropy for the equation of polarized radiation transfer in a homogeneous 
anisotropic Universe was obtained in |lf)j | . The degree of polarization of the background radiation is very sensitive to 
the recombination dynamics and to the reheating epoch. For a Bianchi type-II space-time with a frozen-in magnetic 
field the polarization and anisotropy properties of the cosmic microwave background radiation have been studied in 
11]. Faraday -rotation effects, such as the rotation of the linear polarization plane, are found to be independent of 
the spatial curvature effect and are thus the same as in type-I models. 

Almost immediately after the birth of general relativity (GR), more general geometries, with non-metricity and 
torsion, have been proposed by Weyl and Cartan, in order to incorporate in a geometric framework the effects of 
the electromagnetism and of the angular momentum. Later one, these extensions of the classical GR have been 
incorporated in the different gauge theoretical formulations of gravity |12|. There are different gauge theories of 
gravity in dependence of the choice of the gravitational gauge group and of the gravitational Lagrangian. One of the 
most studied gauge theory of gravity is the Poincare gauge theory, which requires a generalization of the Riemannian 
geometry and the introduction of the torsion of the space-time [Tjj. Non-symmetric gravitation theories, in which the 
gravitational field is described by a non-symmetric metric tensor, whose antisymmetric part couples directly to the 
electromagnetic field, have also been proposed [l4| . 

The study of the propagation of light in a gravitational field has provided the classic observational tests of general 
relativity. The modifications of the geometry due to the non-Riemannian effects, like torsion and non-metricity, 
further affect the motion of photons, by inducing a birefringence of the vacuum which can modify the polarization of 
the photons as they propagate through the gravitational field. Thus the study of the polarization of photons from 
astrophysical or cosmological sources can provide a valuable tool for discriminating between the different modifications 
and extensions of general relativity. Ni 15] has shown that non-metric gravitational fields can single out linear 
polarization states of light that propagate with different speeds. He also proposed to use pulsar polarization data to 
impose constraints on non-metric theories. Limits on the gravity-induced polarization of the Zccman components of 
the solar spectral lines have been obtained in |l6j ] . The polarization of light from magnetic white dwarfs can be used 
to impose constraints on the gravity-induced birefringence of space in nonmetric gravitational theories |17| . By using 
the magnetic white dwarf Grw+70° 8247 polarization data one can obtain the constraint I 2 < (4.9 km) 2 , where /* is 
the non-symmetric gravitational theory charge. 

In conventional Maxwell-Lorentz electrodynamics, the propagation of light is influenced by the metric only, and not 
by the torsion T of the space-time. However, there is a possibility of interaction between light and torsion if the latter 
is non- minimally coupled to the electromagnetic field F by means of a Lagrangian of the form ~ l 2 F 2 T 2 , where I is 
a coupling constant. Several such couplings have been proposed and analyzed recently. A supplementary Lagrangian 
of the form L — XqI 2 (T a A F) T a A F was considered in [TJ|, and it was shown that it can yield birefringence in the 
vacuum. Non-minimally coupled homogeneous and isotropic torsion field in a Friedmann-Robertson- Walker geometry 
affects the speed of light, with the photons propagating with a torsion-dependent speed. In fact, torsion generates 
three major effects affecting the propagation of light: it produces an axion field that induces an optical activity 
into space-time, modifies the light cone structure that yields birefringence of the vacuum and modifies the speed 
of light in a torsion-dependent way [19J. An example of a mctric-affinc gauge theory of gravity in which torsion 
couples nonminimally to the electromagnetic field was considered in [20j. The coupling causes a phase difference to 
accumulate between different polarization states of light as they propagate through the metric-affine gravitational 
field. The model has been constrained by using solar spectropolarimetric observations, which allow to set an upper 
bound on the relevant coupling constant k, k 2 < (2.5km) 2 . 

The confirmation that at least some gamma-ray bursts (GRBs) are indeed at cosmological distances raises the 
possibility that observations of these could provide interesting constraints on the fundamental laws of physics (for 
reviews on GRBs see |2l]-|24|). The fine-scale time structure and hard spectra of GRB emissions are very sensitive to 
the possible dispersion of electromagnetic waves in vacuo. Hence the study of short-duration photon bursts propagating 
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over cosmological distances is the most promising: way to probe the quantum gravitational and/or the effects related 
to the existence of extra-dimensions 25]. The modification of the group velocity of the photons by the quantum 
effects would affect the simultaneity of the arrival times of photons with different energies. Thus, given a distant, 
transient source of photons, one could measure the differences in the arrival times of sharp transitions in the signals 
in different energy bands and use this information to constrain quantum gravity and/or multidimensional effects |26| . 

The announcement of the results of the polarization measurements of the prompt gamma-rays and of the optical 
afterglows of gamma-ray bursts attracted a great interests recently. Coburn and Boggs |2?| claimed the measurement 
of a very large linear polarization, 80 ± 20%, of the prompt gamma-rays from an extremely bright burst, GRB 021206. 
However, this polarization measurement has been criticized by Rutledge and Fox 1281 . who obtained an upper limit 
of only < 4.1% at 90% confidence from the same data, while in Wigger et al. |2!| a value of 41% (+57% -44%) is 
found. One the other hand, the polarization measurements of the optical afterglows are much more convincing. The 
first positive detection of the polarization was for GRB 990510 with the degree of 1.7±0.2%, and since then polarized 
emission has been measured in several other afterglows (see |30] for a review) . 

It is the purpose of the present paper to generalize the transport equations for polarized radiation 0], Q, 
for the case of the polarized radiation propagating in a space-time manifold with non-Riemannian geometry. More 
exactly, we shall consider the case of the so-called Weyl-Cartan s pac e-times, whose geometrical structure is described 
by three tensors: the metric, the torsion and the non-metricity |l3(. In this case the propagation equations for the 
radiation distribution function, for the Stokes parameters, for the linear polarization and for the degree of the total 
polarization can be formulated as a system of a first order partial differential equations on the photon phase-space, 
with the non-Riemannian effects included via the contorsion tensor. 

As an astrophysical application of the obtained equations we consider the propagation of gamma ray bursts in a 
non-Riemannian geometrical background. The general solution of the propagation equation for the total polarization 
degree can be obtained in an exact form. The presence of a non-Riemannian geometrical background, in which the 
electromagnetic fields couple to torsion and/or non-metricity affect the polarization of photon beams. The polarization 
of the beam is generally a function of time, with non-Riemannian effects induced via the contorsion tensor. Since most 
of the gamma ray bursts have a cosmological origin, due to the long propagation times, the influence of the torsion and 
non-metricity could significantly affect their polarization. Therefore we suggest that the observed polarization of the 
prompt gamma ray bursts and of their optical afterglows could also contain a propagation effect component, due to 
a torsion/non-metricity induced birefringence of the vacuum. If such a component could be unambiguously detected, 
this would provide a significant test for the existence of the non-Riemannian geometrical effects in our Universe. 
However, there are many conventional astrophysical mechanisms that could explain the polarization of the gamma 
ray bursts, and affect the propagation of the electromagnetic radiation. The uncertainties in the knowledge of the 
astrophysical environment in which photons propagate, and the difficulty in separating the different physical effects 
contributing to the polarization of the gamma ray bursts make the practical detection of the torsion or non-metricity 
an extremely difficult and very challenging observational task. 

The present paper is organized as follows. The transport equations for polarized light in Weyl-Cartan space-times 
are derived in Section II. The general solution for the total polarization degree for a homogeneous and isotropic 
background geometry is obtained in Section III. Constraints on torsion and non-metricity obtained from the observed 
polarization of gamma ray bursts are obtained in Section IV. In Section V we discuss and conclude our results. 

II. TRANSPORT EQUATIONS FOR POLARIZED LIGHT IN WEYL-CARTAN SPACE-TIMES 

Relativistic transport theory finds its most elegant and natural expression in terms of geometric structures defined 
in the tangent bundle over the space-time manifold. 

Consider a time-oriented Lorentzian four-dimensional space-time manifold M, with metric g of signature 
(+,—,—,—). The tangent bundle T(M) is a real vector bundle whose fibers at a point x £ M is given by the 
tangent space T X (M). If X is a tangent vector and s is a section of T(M), then a connection V is a rule 
Vx (s) for taking the directional derivative of s in the direction X satisfying the properties of linearity in s and 
X, behaving like a first-order differential operator and being tensorial in X. The curvature operator is defined as 
R(X,Y) = V X V Y (s) - VyVx (s) - V[x,Y] (s), where R (^f^, (e*) = ejR^, where {ei,e 2 , ...,e fc } is a local 
frame defined on each neighborhood U C M . The curvature operator satisfies the properties of multi-linearity, anti- 
symmetry and tensoriality. The torsion operator on T(M) is defined as T(X,Y) = VjF — VyX — [X, Y]. This 
is a vector field with components T (t^t, ^t) = (T* u - T^) = T^gfx, where T* v are the Christoffel symbols 
on T(M) and is the torsion tensor |3l|]. Moreover, we assume that the space-time manifold M is equipped with 
a conformal structure, i.e. with a class [g] of conformally equivalent Lorentz metrics (and not a definite metric as 
in general relativity). This corresponds to the requirement that it should only be possible to compare lengths at 
one and the same world point. We also suppose that the connection V respects the conformal structure. Differ- 
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entially this means that for any g G [g] the covariant derivative Vg should be proportional to g: Vg = — 2A (g) g 
\9nv = —ZAxg^ = —Qx^v), where A = A^dx^ is a differential f-form and Q is called the non-metricity [32] • A 
change of calibration of the metric induces a gauge transformation for A: g — > exp (2A) g, A — > A — <iA. Therefore in 
terms of the torsion and non-metricity the connection on M can be expressed as 



r£„ - t£„ + r A „ + A + t„ m A + ~ (q^ a - - q v \) , (1) 

it" + Vi ci f~^. Vi T"i 0^i~»flRol ci rmn r> r~i m t-\i i -fori 4Vr\m +■ Vi cj m nf vi V\tr ncinrr -fV> a <r*2in oro 1 vol ofixriof 1p nracpri r"\+~ i r\n 



where 7^ is the Christoffel symbol computed from the metric g by using the general relativistic prescription, 7^0 

^A A _ a w i, prp Ej- A 

IMV IflU • Ii I/ J U' W11CIC /l^, 



<7 M " {d a g v p + dpg va — d v g a p). Generally, we may represent the connection as = 7* — if A , where K*, the 



contorsion tensor, is a function of the torsion T and of the non-metricity Q |13| 

In the space-time M the instantaneous state of a photon is given by a four-momentum p £ T x (M) at an event 
a; G ill. The one-particle phase space P 7 is a subset of the tangent bundle given by [f| 

P 7 := I x G M,p G (M) ,p 2 = 0} . (2) 

A state of a multiphoton system is described by a continuous, non- negative distribution function f (x,p), defined 
on P 7 , and which gives the number dN of the particles of the system which cross a certain space-like volume dV at x, 
and whose 4- momenta p lie within a corresponding three-surface element dp in the momentum space. The mean value 
of / gives the average number of occupied photon states (x,p). Macroscopic, observable quantities can be defined as 
moments of /. 

Let {x a }, a = 0, 1, 2, 3 be a local coordinate system in M, defined in some open set U C M. Then {d/dx a } is 
the corresponding natural basis for tangent vectors. We express each tangent vector p in U in terms of this basis 
as p = p Q gf^ and define a system of local coordinates {z' 4 }, A = 0, 7 in Tjj{M) as z a = x Q , z a+4 = p a . This 
defines a natural basis in the tangent space given by {d/dz A } — {d/dx a ,d/dp a }. A vertical vector field over TM 
is given by ir = p a d/dp a . The geodesic flow field a, which can be constructed over the tangent bundle, is defined as 
a = p a d/dx a — p a p 1 T^ 1 d/dp^ = p a D a , where D a = d/dx a — pT^d/dp 13 . Physically, u describes the phase flow 
for a stream of particles whose motion through space-time is geodesic. 

Therefore the transport equation for the propagation of a photon beam in a curved arbitrary non-Ricmannian 
space-time is given by 

P a ^-P a ^U£)f = <>- (3) 

By taking explicitly into account the decomposition of the connection in an arbitrary non-Ricmannian space-time 
with torsion and non-metricity, Eq. J3J for the photon distribution function takes the form 

To describe polarized or partially polarized radiation, the distribution function has to be generalized to a distribution 
tensor f a p (x,p) = VJ 2 b _ 1 fab^a where e^f 1 (a = 1, 2) two space-like unit vectors orthogonal to the direction of 
propagation. /„& is a Hermitian and positive matrix, satisfying the conditions Tr(/ a f,(x,p) = f(x,p), det/ a fc (x,p) < 
f(x,p) and f a b = fba- The matrix f a b can be parameterized by the Stokes parameters : P 7 — > [—1, 1], i = 1, 2, 3, 
three measurable quantities, which describe polarization as 33] 

fab = \f{x,p) ( * + 1 2 f + + f 2 ) = \f{x,p) (1+Z-a), (5) 

where a = (oi cr 2 03), with <7j the Pauli matrices, and £ = (£1 £2 £3)- 

In a space-like hypersurface element f2 of P 7 at (x,p) the number of photons linearly polarized along the and 
e^ 2 ) axes are i (1 + ^3) /fi and i (1 — ^3) /O. 5 (1 + ^1) is the result of the measurement of the linear polarization 

in a direction at 7r/4 to the axis. Pj, = +Cl) 1 ^ 2 ^ s the degree of linear polarization. The parameter £2, the 

antisymmetric part of f a b, represents the degree of circular polarization. A measure for the right-hand and left-hand 

1/2 

circular polarization is (1 + £2) /2 and (1 — £ 2 ) /2, respectively. The polarization degree is given by P = f J^ =1 £ 4 2 ^ 

The polarization angle is defined as \& = tan (£1/^3) /2 Q- 

The transport equation satisfied by the distribution tensor L, describing polarized radiation, is 0, 

// ' V Q - ^ /^(x, p ) = 0. (6) 
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In a general non-Riemannian space-time with non- vanishing torsion and non-metricity the equation satisfied by the 
distribution tensor takes the form 

pa ^- pap ^w) ^+W| : -P i/>»-l ,f 'i^ +P a K£h„ + P a Kjf, x = 0. (7) 

As in the Riemannian space-time we may introduce a pseudo-orthonormal basis of vectors e a — e a a d a ( a tetrad), and 
the dual basis 9 a = e"dx a of the one-forms. The components e a a and their reciprocals e"satisfy the relations e°e^ = 5% 
and e^ef = 6^, respectively. The metric tensor can be expressed as g ab = e^e^Tjap, where r) a p is the Minkowski metric 
tensor [3l]]. The object of anholonomity Q ab c = e^e^d^e^ , Q a b c — Q a b d 9dc measures the non-commutativity of the 
tetrad basis. The connection expressed in these anholonomic coordinates is |13|, [3l| 

Tafcc = Tafc.9dc = ~^abc + ^bca ~ ^cab ~ K a b c - (8) 

Any tangent vector p at x can be expressed as p = p a e a . If we set e a — e"e a , then p a = e a a p a and p a = e"p a . In 
terms of the tetrad components the transport equations can be expressed in the forms 

p a DJ = 0,p a D a f bc - p a T d b Jd, - P a T d c Jbd = 0, (9) 

where D a = d a - T b ac p c d/dp b . 

As an application of the tetrad formalism briefly described above we consider the transformation from the coordinate 
basis e a — d a to the basis b a of an arbitrary coordinate system, obtained from a Lorentz boost of the coordinate basis, 
so that b a — L^e/s, e a — L^bp, with the transformation matrices L^ a and satisfying the condition L"L^ = 6%. 

The explicit form of the matrix components is given by L\ — u b , = \a a u a and L@ — a a (6^ + ku^u a ), where 

a a = \/9 aa i ^ = \/~ 9 00 an( l k = (l + u ^/— goo) • The components of the four- momentum with respect to b a are 
p a = LppP . The photon four- momentum is described by its frequency, defined as v — —p ■ u — —p a b a ■ u — p° , and 
two angle for the direction of propagation. Hence, the photon four- momentum can be described in the basis {6 Q }as 
p a = v (1, h). To parameterize the unit three- vector n we introduce spherical coordinates in the momentum space so 
that 



n = fi,n 



1 — " " 2 = y/i - /i 2 cos0, n 3 = \f \ - n 2 sin cj>. (10) 



By assuming that the observer places a two-dimensional screen normal to the propagation direction of the photon 
beam in his rest frame, we can perform a further transformation of the basis so that in the new basis d a — D^bp the 
spatial propagation direction of the photon is now the z-axis. The transformation matrix D@ is given by 



/ 1 \ 

vA-M 2 i>_ 

— /icos0 sin0 y/\ — (j, 2 cos (j> 

\0 — fismcf) — cos0 yl — /.t 2 sin <j> j 



(11) 



The components of the four-momentum in the new basis {rf a }are p a = v (1, 0, 0, 1). With respect to the initial basis 
e Q = 9 Q we have d a = A^ep, where A^ a = L^D^. The inverse transformation is given by e a = A@dp. The components 
of the four-momentum in the two basis are related by the transformations p a = A2f>^ = D^p^ . 

With respect to the basis {d a } the transport equation can be written as 

(Ap^ -jPpP^^j j ab +pap0K ap ^-p a r a Jc b -p a lUac+P a KUcb+r = 0, (12) 

where K a g i are the components of the contorsion tensor with respect to {d a }. The new connection coefficients 7^ 
and K a/3 x c&n be calculated by means of the transformations 

7^ = ^(1^+5,^), (13) 



and 



K af) A = A^Al [k v a K vv " + d v A») , (14) 
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respectively. 
By denoting 



it follows that in a non-Riemannian space-time the radiation distribution function and the Stokes parameters satisfy 
the equations 

Df = 0, D£i — 2 (f 2Q + 7| 3 ) 6 + 2 (i^ + K 2 2 3 ) & = 0, (16) 



£>6 = 0, + 2 (7! + 7l 3 ) 6 ~ 2 (jf*, + Kf 3 ) = 0. (17) 

The evolution equations describing the degrees of the linear polarization and the total polarization P are 

DP L = 0,DP = 0. (18) 

In the limit of the zero contorsion K a/3 Z — > 0, corresponding to the transition to the Riemannian geometry, from 
Eqs. I|16l) - (|18|) we recover the transport equations given in 

III. RADIATION TRANSFER IN ISOTROPIC AND HOMOGENEOUS WEYL-CARTAN SPACE-TIMES 

In order to analyze the influence of the non-Riemannian background on the propagation properties of the electro- 
magnetic radiation we adopt the simplifying assumption of an isotropic, homogeneous and flat Friedmann-Robertson- 
Walker type geometry, that is, we assume that the metric of the Universe is given by 

ds 2 = dt 2 - a 2 {t)(dx 2 + dy 2 + dz 2 ) . (19) 

Generally in the comoving reference frame the torsion tensor is determined by two functions of time T(t) and 
T (t), so that T[j a p = (T(t)/l)5 a p and T a gy = (T (t) /l) £ Q/ 3 7 , where I is a coupling constant and £ a /3y is the totally 
antisymmetric tensor in three dimensions [l9| . By supposing that the theory is invariant under space inversions we 
have T (t) = 0. Therefore the only non- vanishing components of the torsion are T(t) = T\ = T 2 2a = T\ a . The 
non-metricity torsion is defined by three functions of time Qi(t),i = 1,2,3, so that, in an anholonomic basis we have 
Ono = O220 = Q330 = Qi(t), Qooo = and Q on = Q022 = Q033 = Qz{t) • 

By denoting K{i) = — [T(t) + Qz{t)], the transport equation of the total polarization degree P of the gamma-rays 
in a non-Riemannian space-time is given by 



d [id J\ — fi 2 cos 4> d y/l — fi 2 sin 4> d fa \ 9 

vl-- K{t) — 



dt a dx a dy a dz \a ) dv 



P = 0. (20) 



We assume that the photon is emitted at a point Q = (xq, yo, zq) with frequency vq in the direction (/iq, (/>q) and is 
observed at the point O — (x,y,z) with the corresponding four-momentum quantities (is, /j,,<j>). 
Eq. 12Uf) represents a first order partial differential equations with the characteristics given by 

adx ady adz dv , . 

dt= = = ; = t --T, (21) 

which can be immediately integrated to give 



yl — (i 2 cos 4>x — jiy = \f 1 — /Zq cos 4>qXo — (XoVo — constant, (22) 



\f \ — [P- sin (j>x — fxz = J 1 — Hq sin </>oXq — HqZq = constant, (23) 



sin^y — cos <j>z = sin</>oyo — cos^o^o = constant, 



(24) 
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In (av) + J K{t)dt = constant, 



(25) 



x — n 



[ — = constant, y — \J\ — y? cos 4> [ — — constant, z — y/l — y? sin <j) [ — = constant. (26) 
J a J a J a 



Hence the general solution of Eq. (|20|l is given by 



P = P {^s/l — y 2 cos 4>x — fxy, \Jl — p? sin^x — yz, sin^y — cos^z, x — fx J 



dt 

■ 

a 



dt , - . r dt 



y - \j\ - ii 2 cos0 —,z- \J\ - m 2 sin c/> / — ,ln(av) + / K(t)dt . (27) 



a 



For K(t) = 0, Eq. (|25|l gives the usual redshift relation av = a^vo an d Eq. (|27|l describes the propagation of 
radiation in the Ricmannian space-time of general relativity. 

In the following we denote by Pq the value of the polarization corresponding to the propagation of the electromag- 
netic radiation in a curved Riemannian space-time, Pq = P\ K ^ =Q . We also make the simplifying assumption that 
generally the dependence of the polarization on the term In (av) + J K(t)dt is linear, so that 



P = Po 



In (av) 



K(t)dt/) 



(28) 



In order to obtain a quantitative characterization of the non-Riemannian effects on the radiation propagation we 
introduce a parameter (5, describing the variation of the polarization due to the propagation effects, and defined as 

S = ^A- (29) 

In terms of the contorsion tensor, by taking into account that at the receiver a = a rec and by denoting the observed 
frequency of the radiation by v rec , we obtain 



1 



In (a rec v r 



K(t)dt. (30) 



The integration over t can be converted to integration over the redshift z, by using the equality dt = 
(dt/da)(da/dz)dz = —dn (z)dz/(l + z), where dn(z) = (a/a)^ 1 = H^ 1 (z), where H is the Hubble function. We 
assume that generally the Hubble function can be written as 



H (z) - H ^n m , (1 + zf + fi A + Q K fif(z), (31) 

where Hq — 3.24 x 10 -18 ft. s _1 , with 0.5 < h < 1, £! TO ,o ~ 0.3 is the present day matter density parameter, J7a sw 0.7 is 
the dark energy density parameter and £Ik,o is the density parameter formally associated with the non-Riemannian 
effects. / (z) describes the time variation of the contorsion. Qk,o smd f(z) are strongly model dependent quantities 
and generally they depend on the assumed functional form for the contorsion K(t). Therefore in terms of the redshift, 
the parameter 6 can be expressed as 

S( Z) = ** f K ^ dz (32) 



In (a rec v rec ) J z (1 + z ) ^ nmi0 (1 + z f + q a + fl Kfi f(z) ' 
where t H = 1/H = 3.09 x 10 17 /?.- 1 s is the Hubble time. 

IV. CONSTRAINING THE TORSION AND THE NON-METRICITY WITH GAMMA RAY BURSTS 



As an application of the formalism developed in the previous Section we consider the possible effects of a non- 
Riemannian structure of the space-time on the propagation of the gamma-ray bursts. Gamma-ray bursts (GRBs) are 
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cosmic gamma ray emissions with typical fluxes of the order of 10~ 5 to 5 x 10 _4 erg cm -2 with the rise time as low as 
1CP 4 s and the duration of bursts from 10 -2 to 10 3 s. The distribution of the bursts is isotropic and they are believed 
to have a cosmological origin, recent observations suggesting that GRBs might originate at extra-galactic distances 
pi). The large inferred distances imply isotropic energy losses as large as 3 x 10 53 erg for GRB 971214 and 3.4 x 10 54 
erg for GRB 990123 |3. 

The widely accepted interpretation of the phenomenology of 7-ray bursts is that the observable effects are due to 
the dissipation of the kinetic energy of a relativistically expanding fireball, whose primal cause is not yet known |23| . 

The proposed models for the energy source involve merger of binary neutron stars, capture of neutron stars by black 
holes, differentially rotating neutron stars or neutron star-quark star conversion etc. ([^ and references therein). 
However, the most popular model involves the violent formation of an approximately one solar mass black hole, 
surrounded by a similarly massive debris torus. The formation of the black hole and debris torus may take place 
through the coalescence of a compact binary or the collapse of a quickly rotating massive stellar core |23 ■ There are 
still many open problems concerning GRBs, from which the most important is the problem of the source of the large 
energy emission during the bursts. 

The recent observations of the polarization of the GRB's are considered to be of fundamental importance for the 
understanding of the nature and properties of these phenomena. The polarization of the afterglows is well established. 
Typically, the polarization of the afterglows is observed to be at the 1 — 3% level, with constant or smoothly variable 
level and position angle when associated with a relatively smooth light curve. The usual explanation for the polarized 
radiation of afterglows is the synchrotron emission from the fireball when some kind of asymmetry is present. There 
are mainly two kinds of asymmetry considered. One class includes the models assuming that ordered magnetic fields 
play a crucial role. The magnetic fields can be either locally ordered, which corresponds to the magnetic domain 
model |35j| . or even entirely aligned within the ejecta, which is magnetized by the central engine |36j |. In addition, 
small regions in which the magnetic field has some degree of order could be amplified by scintillation [37|], or by 
gravitational micro- lensing (3^ . A different class of models postulates that the fireball is collimated. In this case the 
observer likely sees the fireball off-axis, as the probability of being exactly on-axis is vanishingly small. When the 
line of sight makes an angle with the collimation axis, the asymmetry and hence a net polarization may arise. The 
magnetic fields, compressed in the plane normal to the motion, are postulated to be distributed in a random but 
anisotropic manner Is?!. |39j ]. 

A strong gamma-r ay p olarization may indicate a strongly magnetized central engine, either in pure Poynting- 
flux-dominated form |40| or in conventional hydrodynamical form, but with a globally organized magnetic field 
configuration |4l| . Models involving inverse Compton scattering with offset beaming angles can also give rise to large 
degrees of polarization in gamma-rays |42j . There are also indications that the polarization degree and position angle 
may evolve significantly with time for optical transients |43| . Two-component jet models have also been proposed to 
explain the observed polarization of the afterglows 0] . 

Hence still there arc no definite physical models to consistently predict or explain the observed polarization of the 
gamma ray bursts. 

Since a non-Riemannian geometrical structure of the space-time manifold cannot be excluded a priori, the possible 
effect of a non-zero contorsion must also be taken into account in the description of the propagation of the gamma ray 
bursts. In the following we investigate the effect of different choices of the torsion and non-metricity, corresponding 
to different physical and cosmological models, on the parameter S — (P — Pq)/Pq, describing the deviations in the 
polarization of the gamma-ray bursts due to the propagation effects in a Weyl-Cartan geometry. 



The simplest case that could arise in the analysis of the modifications of the polarization of the gamma-ray bursts 
due to propagation effects is the one corresponding to a constant contorsion, with K(z) — Ka ^constant. Moreover, 
we ignore the contributions of the torsion and non-metricity to the density parameter, by taking Qk.q = 0. Therefore 
we obtain for S the simple expression 



A. 



Models with constant contorsion 




For this model the variation of S as a function of z is represented, for different numerical values of the constant 
contorsion Kg, in Fig. 1. 



9 




FIG. 1: Variation, as a function of the redshift z, of the parameter S (in a logarithmic scale) for gamma-ray bursts propagating 
in a Weyl-Cartan geometry with constant contorsion tensor Ka, for different values of Ka- Ka = 10 _14 s _1 (solid curve), 
Ka = 10~ 15 s _1 (dotted curve), Ka = lCP 16 s -1 (dashed curve) and Ka = 10~ 17 s~ 1 (long dashed curve). For the cosmological 
parameters we have adopted the values a rec = 1, fim.o = 0.3 and Qa = 0.7. The frequency at the detector of the gamma ray 
bursts is assumed to be v rEC = 3 x 10 14 s _1 , corresponding to the optical afterglow emission. 



On the other hand the observational knowledge of 8 would allow to impose some constraints on the value of the 
contorsion tensor, with 

K A < g oM<WV«c) ^ (34) 

F(z) 



where we 



denoted F(z) = J* dz/ (1 + z) Jti^o (1 + zf + ^a- 



B. The Weyssenhoff spin fluid case 



In the framework of the Einstein-Cartan theory, the intrinsic angular momentum (the spin) of a particle is introduced 
via a generalization of the structure of the space-time, by assuming that the affine connection is nonsymmetric. The 
torsion contributes to the energy-momentum of a spin fluid which has the form 

Tj£ = (p + P - 2s2 ) u n u » ~ (p ~ s 2 ) 9iiv, (35) 

where p and p are the energy density and the pressure, respectively, and s 2 = s^s^" is the spin. This energy- 
momentum tensor corresponds to a perfect fluid with spin, called the Weyssenhoff fluid model 01 • The macroscopic 
spin tensor may be expressed in terms of the spin density tensor and the four- velocity of the fluid as = u x s^ v . 
For a pressureless fluid (p = 0) we have p — mn and s — hn/2, respectively, where n is the number of particles with 
mass rn per unit volume and s is the spin density. The effect of the spin is dynamically equivalent to introducing 
into the model some additional non- interacting spin fluid with the energy density p s = pos/o 6 |45|. Although the 
contribution of the spin to the dense matter appears to be negligible small, on a large scale it can produces a 
"centrifugal force" which is able to prevent the occurrence of singularities in cosmology |13|. 

By taking into account that for a Weyssenhoff fluid dominated universe the contorsion is proportional to the spin, 
the parameter S becomes 

, t H K B f° (l + zfdz 

S ( z ) = - uT77. °\ / / (36) 



ln(Wrec)i/z a/Wo (i + zf + n sfi (i + zf + r» A ' 



where Q s .o = p s .o/3Hq is the density parameter of the spin fluid, with p s .o = ft 2 n(0)/I6. Sq can be related to the 
mass to of the particles forming the spinning fluid by means of the relation Kb = p s ,o / v2m. 

The variation of the parameter 5 for a Weyssenhoff fluid filled universe is represented, for different values of the 
cosmological parameters, in Fig. 2. 



The study of the polarization of the gamma ray bursts can impose some constraints on s. Once the functional forms 
of the polarization of the gamma ray bursts at the source and observer are known, from Eq. (|36fl one can evaluate 
the basic physical parameters characterizing the Weyssenhoff fluid and obtain some restrictions on the spin density 
of the particles. 
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FIG. 2: The parameter 5 as a function of the redshift z for gamma-ray bursts propagating in a Weyssenhoff spinning fluid 
filled universe, with the torsion tensor proportional to the spin density, for different values of the cosmological parameters: 
f2m,o = 0.25, £Ia = 0.7 and Q a ,o = 0.05 (solid curve), fi m ,o = 0.20, £Ia = 0.7 and fi s ,o = 0.1 (dotted curve), Q m ,o = 0.15, 
= 0.7 and fi s ,o = 0.15 (short dashed curve) and f2 m ,o = 0.20, Qa = 0.6 and Q s ,o = 0.20 (long dashed curve). The frequency 
at the detector of the gamma ray bursts is assumed to be v rec — 3 x 10 s _1 , corresponding to the optical afterglow emission. 
For the constant Kb we have adopted the value Kb = 2.5 x 10 _16 s _1 



C. Models with modified double duality ansatz 



An extensive study of the spatially homogeneous and S0(3)-isotropic cosmologies in the framework of the 10 
parameter Lagrangian of the Poincare gauge theory was performed in [46(. By evaluating the action in the tangent 
spaces induced by the isotropy group it follows that the torsion tensor has two non-vanishing components h(t) and 
f(t). By excluding the parity violating terms in the Lagrangian the field equations are invariant under the discrete 
transformation / — > — / and s — > — s, where s is the spin scalar. This symmetry rules out classical spin distributions 
like the Weyssenhoff fluid we have previously discussed. It is difficult to find a physical interpretation of a spin tensor 
within this model and therefore the case of a vanishing spin was considered in detail. 

A simple class of solutions of the Einstein-Cartan field equations can be obtained by assuming that the axial 
torsion is non-vanishing and the modified double duality ansatz $ = $o =constant and $ = is satisfied, where 
$ = M + N = R/6 and # = fH + F = R^-yse ^ 6 , with H = h + a/a, M = H + (a/a) H, N = H 2 + k/a 2 - / 2 /4 

/ + (a/a) / . By further imposing the constraints H = and k = 0, it follows that the only non- 



and F 



vanishing component of the torsion is h — —a/a, and the field equations reduce to (a/a) 2 = (47rG/9c4) p — C0/6C4, 
cq,C4 ^constants, and the conservation law for the energy and pressure. In this model the function K(t), describing 
the effect of the torsion on the propagation of the gamma ray bursts is given by K(t) = —Kc^n.a(t), where Kq 
is an integration constant. In terms of the redshift we have K(z) = fsTcln(l + z). The constants —cq/Qc^ can be 
interpreted as an effective cosmological constant A e yj, generated by the presence of the torsion. Therefore in this 
model the parameter 5 (z) is given by 

5(z) = -pp- ^ 1 + *> & (37) 



hx(u rec )j (1 + z) ^/ nmfi (1 + z f + nA 



where Oa = K e ff/p c is the mass density parameter associated to the torsion- generated cosmological constant. 
The variation of 5 as a function of the redshift is presented in Fig. 3. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have generalized the radiation propagation equations for polarized photons in curved 
space-times from Riemannian geometries to Weyl-Cartan type geometries, with torsion and non-metricity. The non- 
Riemannian background enters in the equations via the contorsion tensor and, for particles propagating on cosmological 
distances, can induce a significant and observationally detectable change in the photon polarization state. Of course, to 
fully evaluate this effect one should know the behavior of the contorsion in different classes of generalized gravitation 
theories with torsion. For the Poincare gauge theory and for vanishing non-metricity spatially homogeneous and 
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FIG. 3: The parameter S as a function of the redshift z for gamma-ray bursts propagating in a Riemann-Cartan space- 
time, with vanishing spin density and the modified double duality ansatz satisfied, for different values of the constant Kc- 
K c = 2 x 10~ 15 s _1 (solid curve), K c = lO" 1 ^" 1 (dotted curve), Kc = 8 x 10 _16 s _1 (dashed curve) and K c = 5 x lO^s -1 
(long dashed curve). For the cosmological parameters we have adopted the values a rec — 1, Q m .o = 0.3 and Q,a = 0.7. The 
frequency at the detector of the gamma ray bursts is assumed to be u rec = 3 x 10 14 s _1 , corresponding to the optical afterglow 
emission. 

S'0(3)-isotropic exact solutions for a 10-parameter Lagrangian have been studied in 4jjjj, while the effect of the non- 
metricity on the cosmological evolution has been considered in [84l | . All these solutions are strongly model-dependent 
and generally have a complicated mathematical form. 

On the other hand the study of the polarization of the prompt cosmological gamma ray bursts and of their afterglows 
offers, at least in principle, the possibility of observationally testing Poincare type gauge theories. 

Assuming, as an extreme and perhaps unrealistic case, that the initial gamma ray emission is totally unpolarized, 
as predicted by some models of the gamma ray burst emission (see [2l]|-[24| and references therein), the detected 
polarization amount could be mainly due to the propagation effects of the photons travelling on cosmological distances 
in a non-Riemannian geometrical background. The presence of the contorsion in the radiation transport equations in 
curved space-times could generate a supplementary polarization of the photon beam. If the polarization of the gamma 
ray bursts is indeed due to the coupling between photons and geometry, then the observed degree of the polarization 
must be redshift-dependent, since photons travelling on longer distances will be more affected. The smallness of the 
observed effect even for high redshift sources (P = 1.7 ± 0.2%) also seems to suggest mainly a propagation effect. 

A second important effect, which follows from the consideration of the non-Riemannian structure of the space-time 
is related to the prediction of the frequency dependence of the polarization of gamma ray bursts. Since the dependence 
of the frequency is logarithmic, the total polarization degree of the gamma ray component should be smaller by a factor 
of around 2 than the polarization degree of the optical afterglow emission. Together with the cosmological redshift 
dependence, the frequency dependence of the polarization degree of gamma ray bursts provides a clear observational 
signature of the possible presence of non-Riemannian, Weyl-Cartan type geometrical features in our Universe. 

On the other hand, if the initial polarization degree of the gamma ray bursts could be exactly predicted by some 
photon emission models, then the comparison of the polarization of the gamma ray bursts at the observer and 
emitter could impose some strong constraints on the deviations of the geometry of the Universe with respect to the 
Riemannian background. The torsional effects generated by the spin fluid modify the background geometry and affects 
the propagation of the photons. The study of the polarization of the gamma ray bursts could also impose some limits 
on the spin density and overall rotation of the Universe. 

However, the practical implementation of an observational program aiming at detecting torsion and non-metricity 
from the polarization of the gamma ray bursts could be extremely difficult. From the present observational point of 
view, the polarization measurements of the gamma ray bursts and the interpretation of the data are still controversial. 
There is no general agreement on the detected value of the polarization degree (2?|], l2i |. I 2H . and, since at the moment 
the data quality is insufficient to constrain the polarization degree in most of the cases |29f , a major improvement in the 
observational techniques is required. On the other hand, there are many astrophysical mechanisms which can induce 
an initial polarization of the photons from the cosmological gamma ray bursts, which makes the task of distinguishing 
between the different involved physical processes extremely difficult. The propagation of the photons takes place on 
cosmological distances in a cosmic environment whose properties also present many uncertainties. However, when a 
larger amount of high precision spectropolarimetric data from gamma ray bursts will be available, the possibility of 
testing the foundations of general relativity and other more general gravitational theories could become a reality. 
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